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We present a microscopic theory that shows the importance of spin-orbit coupling in multiferroic 
compounds with heavy ions. In BiFeOs (BFO) the spin-orbit coupling at the bismuth ion sites 
results in a special kind of magnetic anisotropy that is linear in the applied _B-field. We show 
how this interaction is capable of disrupting the magnetic cycloid state of bulk BFO, leading to a 
remarkable level of S-field control of magnetism. 

PACS numbers: 75.85.+t, 71.70.Ej, 75.30.Gw, 77.80.Fm 



The ability to control magnetism using electric fields 
is of great fundamental and practical interest. It may 
allow the development of ideal magnetic memories with 
electric write and magnetic read capabilities [H], as well 
as logic devices based on spin waves that dissipate much 
less energy 0. Usually, the interactions that couple spin 
to electric degrees of freedom are too weak to induce qual- 
itative changes to magnetic states. Hence, the search for 
enhanced couplings has focused on multiferroic materials 
with coexisting magnetic and ferroelectric phases 0. 

The key physical interaction that enables cross cor- 
relation between electric and magnetic degrees of free- 
dom is the linear magnetoelectric effect (LME). It oc- 
curs when the coupling between spin and charge is lin- 
ear in either the external electric field E or the inter- 
nal polarization P, and is quadratic in electron spin. 
In a large class of multiferroic materials, the dominant 
form of LME was found to be due to the spin-current ef- 
fect Q, that couples localized spins according to Hsc = 
J2i< 3 Jsc(P x Rij) ■ (Si x Sj), with Rij the vector link- 
ing the atomic location of spin Si to the atomic location 
of spin Sj. In manganese-based multiferroics, the spin- 
current interaction leads to magnetic induced ferroelec- 
tricity and thus allows magnetic field control of ferroelec- 
tricity 0). 

For -E-field control of magnetism, research has been 
centered instead on iron-based multiferroics, with bis- 
muth ferrite [BiFc03 or BFO] being the most notable 
example Q. At temperatures below 1143 K, BFO de- 
velops a strong electric polarization P = 100 /iC/cm 2 
that points along one of the eight cube diagonals of its 
unit cell [Fig. QJa)]. It becomes an antiferromagnct be- 
low 643 K with Fe spins forming a spiral of the cy- 
cloid type, described by antiferromagnetic Neel vector 
L = sin (q ■ r)q + cos (q ■ r)P. The microscopic origin of 
the cycloid can also be understood as arising from the 
spin-current interaction Q. Plugging the cycloid L into 
JisCi one finds that the lowest energy state is always 
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FIG. 1: (color online) (a) Conventional unit cell for BFO. 
The simple cubic axis x, y, z are denoted by grey vectors. 
The ferroelectric polarization P is shown pointing along [111], 
arising mostly from the displacement of Bi ions with respect 
to the oxygens. The directions X and Y denoted by black 
vectors describe the plane perpendicular to P [T3| ■ (b) Energy 
level diagram for Fe 3+ = [Ar]3d 5 and Bi 3+ = [Pt]6s 2 orbitals 
in BFO; the Fermi level lies just above Eza.- 



achieved when the cycloid wavevector q is perpendicular 
to P. Hence the spins are pinned to the plane formed by 
P and q. This fact has been central to all demonstra- 
tions of S-field control of magnetism published to date; 
the application of an S-ficld poles P from one cube diag- 
onal to another, forcing the spin cycloidal arrangement 
to move into another plane 

A recent experiment suggested that the spin-current 
interaction is not the only LME present in BFO. The 
application of an external -E-field to bulk BFO was shown 
to result in a giant shift of magnon frequencies that was 
linear in E and 10 5 times larger than any other known 
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_E-ficld effect on magnon spectra. 

In this letter, we present a microscopic theory of E- 
field induced magnetic anisotropy, and argue that it can 
provide an effective source of LME in magnetoelectric 
materials with large spin-orbit coupling. Our predicted 
LME explains the origin of the -E-field effect on magnon 
spectra measured in bulk BFO Q. Moreover, we show 
that this effect is capable of qualitatively changing the cy- 
cloidal spin state of BFO, leading to homogeneous mag- 
netism with orientation controllable by the direction of 
the applied i?-ficld. 

Model and microscopic calculation of LME. — Our mi- 
croscopic Hamiltonian for coupling between spin and 
electric-field has three contributions, % = Hiatt + Horb + 
Jiso- The first term is due to the lattice potential, 
Hiatt = -^ 2 V 2 / (2m e ) + V aysta x(r), with m e the free elec- 
tron mass. The second term is called an orbital term 
and is given by H or h = —er-E, with e < and r the 
electron's charge and coordinate, respectively. The last 
and most important term is the spin-orbit interaction, 
"Hso = C,t- • c% with £ = —ir x V the electron's orbital 
angular momentum and <r its spin operator. We take the 
spin-orbit interaction to be dominated by the heaviest ion 
of the lattice, and take bismuth in BFO as a prototypical 
example. 

Single ion anisotropy is known to arise as a correction 
to the total spin energy that is second order in the spin- 
orbit interaction (llj . In our case the largest contribution 
arises in the fourth order of our total Hamiltonian, i.e., 
second order in Hso and second order in Hintt or T-L rb- 
An explicit calculation yields 
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where the spin operator S = J2i=i °i represents all five 
electron spins in the Fe 3+ d-shcll. The numerator of 
Eq. (fTJ) is an outer product between vectors 
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involving 6p and 3d localized orbitals at Bi and Fe, re- 
spectively, with a sum over all vectors R Bl linking the 
central Fe to each of its eight neighboring Bi. We eval- 
uate Eq. ((5]) by taking as Bi orbitals the states \6p x ), 
\6p y ), \6pz), and as Fe orbitals the e g states |3c?3 Z 2_ r 2) 
and |3da;2_ y 2), written with respect to the cubic axis 
x,y,z of BFO's parent perovskite lattice. The Fe t^g 
states need not be considered because they are sa 5 eV 
lower in energy [ill ], and only give a small correction to 
Eq. (TTJ. 

We now turn to an explicit evaluation of the matrix 
elements appearing in Eq. (J5|). The spin-orbit matrix 
element is given by 

(6p„i\(£\6p n ) = -iij h' x n, (3) 
with j] = 0.86 eV chosen to match the spin-orbit split- 
ting measured in isolated Bi ions (l2| . Using symmetry, 



all lattice matrix elements (3d m \ H\ a tt \ 6p n ') can be 
expressed in terms of the direction cosines of i?Bi plus 
only two parameters: V p da = (3d 3z /2_ r 2|Hiatt|6pz') and 
V P dix = (3d x ' z '\H\ a ,tt\6px'), with z' pointing along R Bl . A 
similar procedure can be applied to the orbital matrix el- 
ements (3d m | H or b | 6p„'), reducing them to expressions 
that depend on the direction cosines of i?Bi plus matrix 
elements like (3d x ' z > \z'\6p x >), etc. 

In order to compute the vectors in Eq. ([5]), we need 
to sum over all Bi neighbors forming a distorted cube 
around each Fe. We do this by converting the sum into 
an angular integral, 

I ^latt, orb | 6pn>) « ^ / dfl R [1 + SR ■ V R ] 
X(3d m | Klatt, orb I 6Pn'), (4) 



with 5 R = yR\\P + u±E±j denoting the deviation of 
the Bi ions from the perfect cube. This includes Bi dis- 
placement along P causing ferroelectricity; the displace- 
ment is given by R\\ = 0.116i? B i with R Bi = 4.88 A 0. 
The component of E along P can be neglected (it can 
not compete with the internal field generated by ferro- 
electricity), so we write the external E- field as E± = 
E±[cos (i[j)X + sin(y>)Y"], with the rhombohedral axis 
X and Y shown in Fig. [lja). This perpendicular com- 
ponent induces additional lattice displacement u±Ej_; 
an estimate based on infrared spectroscopy [H| yields 
u ± E ± /R Bi = 2.4 x lCT 4 i; ± /(10 5 V/cm). 

After computing the averages over all matrix elements 
we evaluate the magnetic anisotropy [Eq. ([1])] to give rise 
to two terms: 



H 2 = -~ S-P 
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n E = 



cos {ip)S x + cos ( if} — ^- )S. 



(0) 



Equation (|6]) depends linearly on E±, i.e., it gives rise to 
the LME. 

Even in the absence of an external -E-field, we find a 
magnetic anisotropy, 
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with a coupling energy related to the lack of inversion 
symmetry along P, 
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Taking (E 6p — E 3 d) to be equal to BFO's band gap of 2.8 eV [ijl, and using the tabulated values for the Fe-Bi 
bond V pcl a = -71 mcV and V pdv = -41 mcV [H, we get V\\ = -11 meV and a = 32 fieV= 0.4 K. 

The effect of the external E'-ficld is to introduce magnetoelectric coupling with reduced symmetry; from Eq. @ we 
separate orbital and lattice contributions. The orbital LME is given by 

£orb = (jfrj e^(3d 32 /2„ r 2|z'|6p z /) + V3{3dx'z'\z'\6px') + 3\/3(3dj,/*' WfiPz') + •^={3d X 'y'\y'\6px') 

-(3d 3z ,2^ r 2\x'\6p x ,) + -j=(3d x '2_ y ,2\x'\6p x ')), (9) 



while the lattice LME is 
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Note how these arc physically distinct mechanisms: The 
lattice mechanism is proportional to u±E±, i.e., it arises 
from Enfield induced lattice displacement contained in 
"Hiatt- Plugging the tabulated values for V p da and V p dn 
we get fi att = -5 X 10~ 2 /ieV/(10 5 V/cm). The orbital 
mechanism is instead related to -E-field induced orbital 
admixture, and its matrix elements are not tabulated. 
Assuming (3c£|a^|6p) ~ i?Bi we get an order of magnitude 
estimate of £ orb ~ +30 ^eV/(10 5 V/cm). 

Comparison to experiments. — The experiment in 
Ref. [H discovered a strong dependence of magnon fre- 
quencies on the external -E-field, and used group theory 
to fit two kinds of _E-ficld induced anisotropy: These were 
Fi = r (S/4)(E ± -S)(S-P), andF 2 = -(£/4)E ± -[(S$- 
Sx)X + (2SxSy)Y}. It was shown that only F 2 would 
give rise to the observed linear in E± magnon shift, and a 
fit of £cx P = +55 ^cV/(10 5 V/cm) with a = was estab- 
lished at T = 300 K. To compare this result to our theory, 
we write our Eq. ([6]) in the rhombohedral basis and get 
that it is equal to F 2 + 2\/2i r i- Thus our Eq. © can be 
expressed as a function of the two anisotropy terms of 
Ref. and explains the origin of the interaction leading 
to electrical control of magnons in BFO. 

Our calculated zero-field anisotropy energy a = 32 /xeV 
is close to the value of a « 20 fieV extracted from neutron 
diffraction experiments (l7| . 

Electric-field control of magnetism. — To find out 
whether our effect can be used to control magnetism us- 
ing an external E- field, we incorporate Eq. (pt into the 
usual continuum free energy model for BFO [18l420| . 
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FIG. 2: Electric-field induced magnetic phase diagram for 
BFO. Here Ex and Ey are projections of the external E- 
field in the plane perpendicular to P; the axis X, Y are 
shown in Fig. [TJa). A transition from cycloid to homogeneous 
magnetism is predicted for certain directions in the X — Y 
plane. In this case the Neel vector L will point along one 
of the conventional cubic directions h = x,y,z, and due to 
the Dzyaloshinskii-Moriya coupling the magnetization M will 
point along P x L. Thus an external _B-field is able to control 
the direction of the two magnetic order parameters M and 



Here L is the Neel vector, and the first and second terms 
inside the brackets of Eq. (fTTj) arise from the exchange 
interaction between spins; the third term arises from the 
continuum limit of the spin-current coupling, leading to 
a' = J sc (Sl /2) 5 / 3 /(2S'/iB) 2 , with Q = 124.32 A 3 the 
unit cell volume in BFO. This term is often called the 
flexoelectric interaction, and explains the origin of the 
cycloid in BFO when c'q 2 w 1 18 
term is the continuum limit of Eq 

(n £)/(2Sv B ) 2 . 
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The minimum free energy state L(r) can be found us- 
ing functional derivatives in the same way as done in 
Refs. 18, [l9[. The result is summarized in Fig. [2] As 
the electric field is increased, the anisotropy energy fa- 
vors an anharmonic cycloid ground state with L forming 
a square wave along one of the three cubic directions x, 
y, or z, depending on the direction of E±. When E± be- 
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comes larger than a certain critical value, we get a phase 
transition to a homogeneous L, effectively destroying the 
cycloid state. The origin of this phase transition is the 
competition between .E-field induced anisotropy and the 
nexoelectric interaction. The free energy of the cycloid 
state is J" C y C ioid ~ -(£'£a/2)(cos 2 (q • r)) - c'q 2 /2 = 
—£'E±/4: — c'q 2 /2. Compare this to the free energy 
of the homogeneous state, -7-Hom ~ —£'Ej_/2; as E± 
increases, eventually we will have £,'E± > 2c'q 2 and 
Tkom < -^cycloid, inducing a transition to the homoge- 
neous state. A detailed calculation yields the following 
critical field: 



(^X) = T 
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for < ip < 120°, and similar expressions for ip — > 
i> - 120° and ip ip - 240°. Remarkably, the critical 
field is infinite when E points antiparallcl to one of the 
cubic directions. This is also easy to understand from the 
symmetry of Eq. ([6]): for example, when E± points along 
the projection of —x in the X — Y plane (ip = 0°), the 
electric field anisotropy energy is the same for L along y 
or z; thus when L is a cycloid in the y — z plane, it is able 
to simultaneously minimize both the i5-field anisotropy 
and the nexoelectric energies; in this situation it is ener- 
getically favorable for L to remain a cycloid. A similar 
situation applies for E \\ — y or E \\ —z. 

An important point is that P can be poled by the ex- 
ternal -E-field, changing the effective direction of Ej_ in 
Fig. [2] (note that Fig. [5] assumes P || [111] at all magni- 
tudes of E± ) . To avoid poling, one can apply the i?-field 
with a component along the [111] direction. For example, 
using E = £[cos(30°)Z - sin(30°)X] allows control of 
magnetism without changing P, at the expense of having 
E± = E/2. Using £ cxp = 55 ^eV/(10 5 V/cm) we get that 
a minimum E = 1.3 x 10 5 V/cm is required to induce 
the homogeneous state, a value well into the practical 
range. To confirm our theory we propose the application 
of the -E-field to bulk BFO along this specific direction. 
The homogeneous state has as its optical signature the 
presence of only two magnon Raman modes [2l| (the sig- 
nature of the simple antiferromagnet) instead of five or 
more cyclonic magnons [53 ] . 

The ability to switch from cycloidal to homogeneous 
spin order by way of an external i?-field represents a 
significant step towards the development of spin-based 
technologies. It enables the engineering of logic gates 
with magnetic excitations. In addition, BFO is known 
to have a weak magnetization M oc Z x L that is gen- 
erated by the Dzyaloshinskii-Moriya interaction [23. 24 [. 
of magnitude |M| = 0.06/is/Fe [25j. Since the direction 
of M is tied to L, our mechanism allows -E-field control 
of M, a functionality that will be useful for the design 
of magnetic memories that dissipate less energy [l[. 



Conclusions. — We presented a microscopic theory of 
E-Rcld induced magnetic anisotropy, and showed how it 
gives rise to an additional linear magnetoelectric effect 
(LME) in multiferroic materials. The origin of this spe- 
cial kind of LME is based on the combination of two 
factors: The presence of a non-magnetic ion with large 
spin-orbit coupling, and a significant amount of inver- 
sion asymmetry (induced e.g. by ferroelectricity). For 
BFO, the presence of this additional LME implies that 
its magnetic cycloid can be converted into a homoge- 
neous state under the application of a practical external 
i?-ficld. Thus, it shows that _E-field control of magnetism 
in BFO can be much more complete than what has been 
demonstrated so far. 

Our research was supported by the NSERC Discovery 
program. 
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